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We report x-ray reflectivity (XR) and small angle off-specular diffuse scattering (DS) measure- 
ments from the surface of liquid Indium close to its melting point of 156° C. From the XR measure- 
ments we extract the surface structure factor convolved with fluctuations in the height of the liquid 
surface. We present a model to describe DS that takes into account the surface structure factor, 
thermally excited capillary waves and the experimental resolution. The experimentally determined 
DS follows this model with no adjustable parameters, allowing the surface structure factor to be 
deconvolved from the thermally excited height fluctuations. The resulting local electron density 
profile displays exponentially decaying surface induced layering similar to that previously reported 
for Ga and Hg. We compare the details of the local electron density profiles of liquid In, which is a 
nearly free electron metal, and liquid Ga, which is considerably more covalent and shows directional 
bonding in the melt. The oscillatory density profiles have comparable amplitudes in both metals, 
but surface layering decays over a length scale of 3.5 ± 0.6 A for In and 5.5 ± 0.4 A for Ga. Upon 
controlled exposure to oxygen, no oxide monolayer is formed on the liquid In surface, unlike the 
passivating film formed on liquid Gallium. 

PACS numbers: 61.25.Mv 68.10.-m 61.10.-i 
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I. INTRODUCTION 

The structure of the free surface of a liquid metal is 
fundamentally different from that of a dielectric liquid 
such as water or Argon. For dielectric liquids, the inter- 
atomic or intermolecular potential is long-ranged and the 
nature of the interactions does not change significantly 
across the liquid-vapor interface. Theoretical modeling 
and molecular dynamics studies result in an interfacial 
density profile that varies monotonically from the bulk 
liquid density to the bulk vapor densityi This can be 
contrasted with results of computer simulations of dielec- 
tric liquid noble gases in contact with their solid phasen^i 
There, the hard wall provided by the solid phase is shown 
to induce surface layering in the liquid with a well defined 
lamellar structure. 

In liquid metals, the potential energy function changes 
drastically from a screened short-ranged Coulomb po- 
tential in the bulk liquid phase where the conduc- 
tion electrons are delocalized to a long-ranged van 
der Waals type potential in the vapor phase where all 
electrons are localizedi^ In this case, density functional 
calculations^ and computer simulations^ indicate an os- 
cillatory surface-normal density profile. This surface in- 
duced layering can be attributed to the fact that strong 
Coulomb interactions between the abruptly truncated 
conduction electrons at the surface and the ion cores im- 
pose a constraint ordering the surface in a well defined 
lamellar structure. This constraint can be thought of as 
an effective hard wall ionic potential at the surface where 
the conduction electron density abruptly terminates. 

The existence of surface layering in liquid metals has 



only recently been verified unambiguously by experiment, 
first for liquid Hg^ and then for liquid Ga^ Comparison 
of the temperature dependence of this layering, however, 
reveals qualitative differences between the two metalsi^jiS 
This observation indicates that details of the interactions 
such as the degree of covalency in the metal may affect 
the surface structure. 

For crystalline metals that can be described by weak 
pseudo-potentials, electron-ion scattering is relatively 
weak and the non localized or itinerant electrons can be 
treated as if they were essentially freeiii Alkali metals 
with their nearly parabolic energy bands are the prime 
example of this class of nearly free electron (NFE) like 
metals. In analogy with this, the term NFE metals ap- 
plies to liquid metals whose conduction electrons have 
similar itinerant character and are only weakly pertur- 
bated by a small ionic pseudopotential.^^ NFE liquid 
metals can be contrasted with those whose conduction 
electrons are partially localized in covalent bondsii^ Ev- 
idence for such covalent character in liquid metals may 
be obtained from measurements of the optical properties 
as well as from the bulk liquid structure factor^ 

In fact, both liquid Ga and Hg show substantial devi- 
ations from NFE behavior. In liquid Ga, a shoulder in 
the bulk liquid structure factor indicates the presence of 
directional bonding in the meltii^*iS4i2ii& Liquid Hg dis- 
plays a pronounced asymmetry in the first peak of the 
bulk liquid structure factor, and the optical constants of 
liquid Hg deviate markedly from the predictions of the 
Drudc theory for free electronSfi^ which, for example, 
adequately describes the optical constants of the alkali 
metals>22i2i In order to understand which aspects of sur- 



face layering may be affected by deviations from NFE 
character, it is necessary to investigate liquid metals with 
less tendency towards covalent bonding. The alkali met- 
als would be ideal candidates in this regard. Unfortu- 
nately, the very low surface tension of liquid alkali metals 
will result in a large thermally induced surface roughness, 
making it extremely difficult to observe the surface layer- 
ing experimentally. This is exacerbated by the fact that 
alkali metals are very reactive in air but cannot be in- 
vestigated under Ultra High Vacuum (UHV) conditions 
due to their high vapor pressure. Liquid Indium, which 
has a high surface tension, reasonably low melting point 
and low vapor pressure does not suffer from these short- 
comings. Moreover, it is considered to be NFE like for 
the following reasons. First, the first peak in the bulk 
structure factor of liquid In is highly symmetric, indi- 
cating the absence of significant orientational bonding in 
the melt.^^ Second, quantitative analysis of the shape of 
the bulk liquid structure factor yields a twelve-fold co- 
ordination with an interatomic distance that decreases 
with increasing temperaturei^ This is precisely what is 
expected for a liquid that can be described by an ideal 
hard sphere model. Third, the measured optical proper- 
ties of liquid In agree with the Drude free electron theory 
over a large range of wavevectorsi^ Even more interest- 
ing, the optical constants of mixtures of liquid Hg and 
In differ appreciably from free electron behavior for Hg 
rich compositions, but gradually approach free electron 
behavior with increasing In content. ^^ A main objective 
of this study is to compare the surface structure of the 
more NFE like liquid metal In to that of the more cova- 
lent liquid Ga. 

This fundamental interest in the relationship between 
the surface structure and electronic properties of liquid 
metals bears directly on the topic of surface reactions, 
which is of considerable practical importance. This is 
especially true for liquids, where the atoms are mobile, 
and reactions at the surface can directly affect alloying, 
phase formation and other properties of the bulk. For 
this reason, surfaces and reactions at surfaces play a 
critical role in process and extractive metallurgy«2& We 
previously demonstrated that on exposure to oxygen the 
surface of liquid Ga becomes coated with a uniform 5 A 
thick passivating oxide film that protects the underlying 
bulk phase from further oxidationiSi We will show below 
that the oxidation behavior of liquid In is fundamentally 
different from that of Ga. 



II. EXPERIMENTAL DETAILS 

The experiments reported here were primarily per- 
formed at the beamline X25 at the National Synchrotron 
Light Source (NSLS). Following a toroidal focussing mir- 
ror, a Ge (220) crystal was used to select the x-ray wave- 
length A = 0.653 A and to tilt the beam downward onto 
the horizontal liquid metal surface. The geometry of 
this liquid spectrometer has been described elsewhere^ 



Additional data were acquired at beamline X22B at the 
NSLS, using a similar scattering geometry and an x-ray 
wavelength of 1.24 A. The size of the x-ray beam striking 
the liquid metal surface was determined by 1 mm (hor- 
izontal) X 0.1 mm (vertical) slits downstream of the 
monochromator and the measured signal was normalized 
to the incoming beam intensity. 

The experimental resolution, determined primarily by 
slits in front of the detector, was varied for different ex- 
perimental configurations. Specular reflectivity was pri- 
marily measured using a conventional Nal scintillation 
detector. Typical slit settings for measuring this spec- 
ular reflectivity at X25 were 4 mm (h) x 4 mm (v). 
At 600 mm from the sample this yields a q^ resolution 
of 0.06 A~^. For comparable resolution at the longer 
wavelength at X22B, a vertical detector slit of 8 mm was 
used. Small angle off-specular surface diffuse scattering 
was measured using a position sensitive detector (PSD) 
aligned within the reflection planed For a given inci- 
dent angle, the PSD detects both the specular reflection 
and off-specular diffuse intensity in the reflection plane 
up to « 0.5° away from the specular peak. We used a 
coarse vertical resolution for the PSD in the multi chan- 
nel mode (64 channels on 20 mm) . Away from the specu- 
lar peak, where the intensity varies slowly with position, 
data from neighboring channels were averaged together 
for improved statistics. Counting all photons hitting the 
detector, and disregarding the positional information, the 
PSD can be used as a conventional counter. In this "sin- 
gle channel" mode the slit settings were 20 mm vertically 
and 3 mm horizontally, corresponding to a Qz resolution 
of 0.32 A"-'^. For all experimental arrangements, the sig- 
nal originating at the surface was separated from the 
bulk diffuse background by subtracting intensities mea- 
sured with the detector moved approximately ±0.4° (cor- 
responding to « 0.67 A ) transverse to the reflection 
planei2& 

The linearity of the PSD was checked by comparison 
of diffuse scans taken with the PSD to diffuse data taken 
by moving the scintillation detector within the reflection 
plane, and by measuring the specular reflectivity with 
the stationary PSD in a single-channel mode. As we will 
discuss below, once scaled for the different resolution, the 
different types of scans perfectly overlap each other. 

Indium of 99.9995% purity was contained in a Molyb- 
denum pan within an Ultra High Vacuum (UHV) cham- 
ber to maintain a clean surfaced Details of the UHV 
chamber design can be found elsewhereiSS The sam- 
ple was heated above the melting point of liquid In to 
170 ± 5 °C. After a three day bakeout, the pressure was 
8x 10^^° Torr, with a partial oxygen pressure < 8x 10^^^, 
for which the formation time for an oxide monolayer is 
calculated to be several days Ji Although patches of ox- 
ide were present on the In ingot when it was introduced 
into the chamber, this oxide was removed by Argon ion 
sputtering (2.5 keV, 20 /xA ion current) at an Ar par- 
tial pressure of about 5 x 10~^ Torr. The impact of 
the ion beam induced mass flow along the surface that 
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FIG. 1: Sketch of the geometry of x-ray scattering from the 
hquid surface with a and jS denoting incoming and outgoing 
angle, the incoming and outgoing wavevector kin and kout re- 
spectively and the azimuthal angle 29. The momentum trans- 
fer q has an in-plane component qxy and a surface-normal 
component q^. 



transported the oxide patches towards the ion beam, thus 
cleaning the entire liquid-vacuum interface and not only 
the area hit directly by the Ar+ ionsi^ Data presented 
here were always taken within one to two hours of sputter 
cleaning the sample. Since data acquired without sput- 
tering but within a few hours reproduced each other, and 
because of the unique growth properties of oxide on liq- 
uid In which will be presented below, we are confident 
that the measured scattering is intrinsic to a clean liquid 
Indium surface. 

A major improvement over experiments described pre- 
viously is the incorporation of active vibration isolation. 
Previous measurements eliminated mechanically induced 
surface waves by utilizing the viscous drag at the bot- 
tom of thin samples. These samples had curved sur- 
faces, requiring a time consuming technique to measure 
the reflectivitjiSS and making it very difficult to measure 
at grazing incidence, important for the detection of in- 
plane ordering. In this work, the rigid UHV chamber and 
ion pump assembly were mounted on an active vibration 
isolation unit^^ similar to that used previously in our 
studies of liquid Hg.^ With this system we were able to 
use a sample pan 5 mm deep and 60 mm in diameter, re- 
sulting in flat samples where the angular deviation of the 
surface normal from the vertical across the surface was 
negligible compared to the critical angle for total external 
reflection of x-rays. 



III. THEORY 

The scattering geometry is given in Figure ^ defining 
the liquid surface to be lying in the x-y plane, with x- 
rays incident at an angle a and collected by the detector 
at an elevation angle /9 and azimuthal angle 29. The 
momentum transfer q can be decomposed into surface 
normal, q^, and surface parallel, qxy, components given 



q^ = ^ (sin f3 + sin a) 



and 
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(1) 



- ^cos^ a + cos^ /? — 2 COS a cos (3 cos 26 . 



In the following, we develop a formula for the differen- 
tial cross section da/dfl for x-ray scattering from a rough 
and structured liquid surface. The general expression for 
the differential cross section for x-ray scattering from a 
three dimensional electron distribution p(r) at a given lo- 
cation r can be expressed in terms of the electron density 
correlation function {p{f)p{ri)) = {p{r — rf)p{0))i^ 



m=y{^yidHr-ff){pif-ff)x 
xp(0)) exp [iq- (r — r/)] 



(2) 



where V is the illuminated volume. If the density dis- 
tribution is homogeneous within the x-y plane but inho- 
mogeneous normal to the surface, the density correlation 
function depends on the different positions f^y — r^yi on 
the surface and the distances z and zi from the surface. 
For an x-ray beam of a cross sectional area ^g incident 
on the surface at an angle a relative to the x-y plane, the 
illmninated area is v4o/sin(a) and the differential cross 
section can be writteni^ 
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xp(0, zl)) exp [iq^{z - zl) + iq^y ■ r. 
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For a homogeneous liquid that is sufficiently far away 
from any critical region, the only long wavelength excita- 
tions that give rise to significant scattering are thermally 
excited surface capillary waves. In this case, there is some 
length ^ such that for \rxy\ > £, the density density cor- 
relation function can be expressed in terms of a density 
profile p[z — h{fxy)] defined relative to the local position 
of the liquid surface h{fxy) and capillary wave induced 
relative variations of the height, h{fxy) — h{fxyf)^ 

{p{rxy,z)p{0, Zl)) ^'^^^ p{z - [h{fxy) - hm)~p{z') (4) 

Defining 5p{rxy,z) = p{fxy,z) - p{z - h{rxy)), the dif- 
ferential cross section from a liquid surface can now be 
written: 

7S = -4^f-^l [dzdz/d^rxyX 

X exp [iqz{z - zl) + iq^y ■ r^y] x 

X {p{z - [h{rxy) - h{Q)])p{zi) + {5p{rxy, z)5p{Q, zl))} 

(5) 
The second term in the integrand is only non-vanishing 
for regions \rxy\ < C- This term gives rise to both the 



bulk diffuse scattering and surface scattering at momen- 
tum transfers having large qxy coniponentsi^ The effect 
of this term can be separated experimentally from the 
specular reflection and may be omitted in the following 
discussion of scattering close to the specular condition. 
A change of variables yields 
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where Qc is the critical angle for total external reflection 
of x-rays and the surface structure factor 



$(g,) = ^/dzi^exp(zg,z) = 
= 7^ / dzpiz) exp(ig2z) 



(7) 



is the Fourier transform of the local or intrinsic density 
profile p{z), independent of fxy- In this equation, poo is 
the bulk electron density. 

For the case of thermally excited capillary waves on 
a liquid surface, the height fluctuations can be charac- 
terized by their statistical average {\h{'rxy)\'^) -^^ Sinha et 
al. have shown that integration over these height fluctu- 
ations yields the following dependence of the scattering 
on Qxy and q^i^ 
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(i^rry(exp {iqzh{rxy)}) exp [iq^y ■ r^y 



with C to be determined and 
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The value of C can be determined by integration over all 
surface capillary modes having wavevectors smaller than 
the upper wavevector cutoff qmax = tt/^: 



X (exp {iqzh{fxy)}) exp [iq^y ■ Vxy] ~ 
« An'^limf^ ^o{ejip{iqzh{fxy)}) = Aw'^ 



(9) 



Therefore, the properly normalized differential cross sec- 
tion for scattering of x-rays from a liquid surface is: 



da _ j4o /'qc\'^ ksT 
dfl sin ( 
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The intensity measured at a specific scattering vector q 
is obtained by integrating Ea.llUlover the solid angle dfi 
defined by the detector acceptance of the experiment. For 
the specular reflection qxy = (a = /3 and 20 = 0), the 
integral is centered at qxy = 0. The projection of the de- 
tector resolution onto the x-y plane is rectangular'^^ and 



the above mentioned integration has to be done numer- 
ically. However, an analytical formula can be obtained 
for the specular reflectivity if the detector resolution is 
assumed to be a circle of radius g^es , independent of q^^: 
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|^exp[-(T2^gf] 



(11) 



jor qz ^ bqc 



is the Fresnel reflectivity from classical optics for a 
flat {h{fxy) = h{0) for all h{r'xy)) and structureless 
(|$((7z)P = 1) surface. The quantity CTcw denotes the 
capillary wave roughness: 
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In the experiments presented here, the measured specular 
reflectivity was analyzed by numerical integration of the 
theoretical model over the rectangular resolution deflned 
by the detector slits. In practice, the specular reflection 
is measured with coarse resolution and the results based 
on Eq. ^2 fall within 5% of the exact values obtained 
from integrating over Eq. ^| On the other hand, the off 
specular diffuse scattering is measured with much flner 
resolution and near the specular peak, qxy = 0, it has 
to be analyzed by numerical integration of Eq. IIUI over 
the slit defined resolution. In addition, numerical inte- 
gration of Eq. ^1 is necessary to obtain the correct ratio 
of specular reflectivity to off-specular diffuse scattering. 
In Eq.(??) we have introduced a local structure factor 
which is the Fourier transform of the local or intrinsic 
density profile. An alternative representation is to define 
a macroscopically averaged density profile {p{z))t such 
that 



1 

Poc 



dz ; evi-pliqzz) 

dz 



<i>(g,)exp(-l/2aL9.') 



This macroscopically averaged density profile, {p{z))x, 
can be shown to be equal to the convolution of the local 
density profile, p{z), with the associated Gaussian distri- 
bution of height fluctuations characterized by 
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Unfortunately, it is difficult to extract a precisely defined 
{p(z))t from a single reflectivity measurement because 
of the above mentioned problem of the variation of the 
resolution with q^ . 

A second difficulty with reflectivity measurements, 
even when [^(g.)^ can be separated from the thermal 
effects, is that the direct inversion of |$(g2)P in order to 
obtain p{z) is not possible because the phase information 



necessary to perform the Fourier transform in Eq. ?? is 
lost in any intensity measurement. The common practice 
to overcome this problem is to assume a model for the 
density profile, which is then Fourier transformed and 
fitted to the experimentally determined structure factor. 
In this study, the liquid metal is modelled by layers of 
different electron densities parallel to the surface. The 
layers represent planes of atoms, separated by a distance 
d, and having increasing width the further they are below 
the surface. Mathematically, the electron density is de- 
scribed by a semi-infinite sum of Gaussians, normalized 
to the bulk density poo'- 



P{z) = Pc 



E 

n=0 



V2^ 



exp [-{z - ndfl2al] (g) Fi,,{z). 



(13) 

Here (g) denotes convolution, Fin{z) is the atomic x-ray 
scattering form factor of In, and tr^j — na^ + CTq, where 
a and ctq are constants. This form for (T„ produces a 
quadratic increase in the Gaussian width with distance 
z below the surface, so that the parameter a is related 
to the decay length for surface layering and the model 
density approaches the bulk density poc for (t„ 3> d. An 
advantage of using Gaussian functions to model p{z) is 
that an expression for {p{z))t is obtained from Eq. 1131 
by replacing an with a^ and defining (cr^)^ = a-^,^ + a^- 
Note that the effective width of the individual layers be- 
comes (cr^)^ which explicitly demonstrates the additiv- 
ity of the different contributions to the surface roughness 
from capillary modes (for which q^y < Qmax) and short 
wavelength modes (for which q^y > Qmax) that are incor- 
porated into the model for p{z). 



IV. EXPERIMENTAL RESULTS 

The specular refiectivity measured from liquid In at 
170± 5 °C demonstrates the presence of surface layering. 
Figure|2]shows the refiectivity as a function of q^ , normal- 
ized to the Fresnel refiectivity of a sharply terminated In 
surface with a critical angle qc = 0.052 A~^. The circles 
show the data taken with the scintillation detector, with 
a qz resolution of 0.06 A~^. Open circles represent the 
small angle data taken at X22B (A=1.24A) whereas the 
closed circle data have been taken at X25 (A=0.65 A). 

Filled squares represent data taken with the PSD in a 
single channel mode, with the vertical resolution coarser 
by a factor of 5. The large increase in signal with decreas- 
ing resolution immediately demonstrates the significance 
of the diffuse scattering. The reflectivity at small and in- 
termediate wavevectors falls below the Fresnel value cal- 
culated for an ideally flat and abruptly terminated sur- 
face. This is due to the roughness of the liquid In surface 
which scatters photons out of the specular condition. The 
most prominent feature in these data is the broad peak 
centered near 2.2 A"-'^. This peak is due to constructive 
interference of x-rays from layers ordered parallel to the 
surface. This very distinctive structural feature of liquid 
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FIG. 2; Specular x-ray reflectivity for liquid Indium at 170 °C 
taken witli vertical detector resolutions of 0.06A~^ (scintilla- 
tion detector; open circles: X22B, closed circles: X25) and 
0.32A^^ (PSD in single channel mode; squares: X25). The 
reflectivity is normalized to the Fresnel reflectivity Rf of a 
flat In surface. Solid lines are fits as described in the text. 
Data from X22B was not incorporated in any fits. 



metal surfaces has been observed earlier for Hg^ and Gan^ 
The two data sets show the same features, although as 
a consequence of integrating the diffuse scattering over 
coarser resolution, the intensity is greater for data taken 
with the PSD. 

The lines in Figure |21 illustrate fits of the density pro- 
file model discussed above to these data. Three sur- 
face layering parameters, d, a, and ao^ determine the 
form of $(92), which is then used in Eq. 1101 along with 
T = 170°C and the experimentally determined value 
for the surface tension of liquid In at that temperature, 
7 = 0.556 N/m."^^ As discussed above, we assign the short 
wavelength cutoff for capillary waves, qmax ~ t^/S. , to 
be ~ lA where £, was taken to be the nearest neigh- 
bor atomic distance in the bulk mclt.'^^ The integral over 
Eq. 1101 is then performed numerically over the known 
resolution volume. Surface layering parameters obtained 
from the scintillation detector data are d = 2.69±0.05 A, 
a = 0.54 ± 0.06 A, and ctq = 0.35 ± 0.04 A (solid line in 
Figure |2Jl . Essentially the same parameters (with larger 
error due to poorer statistics) result from a fit of the 
model to the PSD data (broken line in Figure |21). 

Although the resolution dependence of the refiectivity 
of liquid metals provides a test of the thermal capillary 
wave theory prediction, a more rigorous test is the mea- 
surement of the spectral density of the off specular dif- 
fuse scattering. We have measured diffuse intensity over 
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FIG. 3: Diffuse scattering as a function of scattering angle 
/3 for different fixed incoming angles a. Solid lines: Dif- 
fuse scattering calculated from the experimentally determined 
structure factor with no further adjustable parameters. In- 
set: linear plot emphasizing the fit near the specular peak for 
a = 4.5° 



a (3 range straddling the specular condition. Intensities 
normalized to the direct beam are shown in Figure 13 for 
several choices of a. Data at extreme values of (/3— a) are 
limited by the intense background due to bulk scattering 
at large /? and by the scattering geometry at small (3. 
The asymmetry of the wings centered around the specu- 
lar ridge {/3—a — 0) arises from the /3 dependence of both 
the exponent j] and the surface structure factor <&((Zz). 

The solid lines are obtained by calculating the inten- 
sity from Eq. IIOI integrating numerically over the reso- 
lution function and subtracting a similar calculation for 
the background scan taken out of the reflection plane. 
The calculation incorporates the surface structure fac- 
tor ^{qz) determined from the reflectivity measurements 
along with the other fixed quantities (7, fc^T, and gmax), 
without further adjustable parameters. We find excel- 
lent agreement between this model and the experimental 
data for the entire range studied, including the wings 
(Figure |2Jl and the specular region (inset of Figure |2J|. 
The implications of this agreement will be discussed in 
the next section. 

To study the oxidation properties of liquid Indium we 
exposed the In surface to controlled amounts of oxygen 



through a bakeable UHV leak valve. During oxidation, 
macroscopic oxide clumps, large enough to be observed 
by eye, formed at the edges of the sample, while the cen- 
ter remained clean. X-ray reflectivity measured during 
oxidation showed no changes until the floating oxide re- 
gions grew large enough to reach the area illuminated 
by the x-ray beam, at which point the macroscopically 
rough surface scattered the reflected signal away from the 
specular condition. This result is in sharp contrast to the 
formation of the highly uniform, 5 A thick, passivating 
oxide layer we previously observed to form on liquid Ga 
under the influence of the same amount of oxygen^ 



V. DISCUSSION 

The macroscopic density profile extracted from the re- 
flectivity measurements is resolution dependent since the 
density oscillations are smeared out by thermally acti- 
vated capillary waves as described in the Theory section. 
This smearing depends on the temperature and surface 
tension and hence also varies for different liquids. The 
macroscopic density profile directly extracted from the 
experiment is the averaged density profile {p{z))t shown 
in previous publicationsiiSi2iS& In order to compare the 
intrinsic layering properties in different liquid metals, it 
is necessary to remove these thermal effects that vary 
from metal to metal and experiment to experiment and 
to obtain the intrinsic or local density profile p{z). Previ- 
ous temperature dependent reflectivity measurements of 
liquid Ga around the specular position have shown that 
the resolution dependent roughness is well described by 
the capillary wave (CW) predictionil. Similar conclusions 
were reached from T-dependent measurements of liquid 
paraffins.^" Here, we employed a different approach to 
verify that the CW prediction holds for liquid In, per- 
forming diffuse scattering measurements at a single tem- 
perature at different angles of incidence. The magnitude 
and the angular dependence of the diffuse scattering is in 
full compliance with the predictions of the capillary wave 
model without using any adjustable parameters. This 
is demonstrated by the perfect agreement of the diffuse 
scattering profiles with the theoretical curves in Figure^ 
This allows us to directly calculate the local density pro- 
file p{z) by simply setting 77 = 0. A comparison between 
the local density profile p{z) and the temperature aver- 
aged density profile {p{z))t is made in FigureEJ^a). 

In this single-T approach the effects of the intrinsic, 
T = 0, roughness ctq and the cutoff g,nax cannot be 
unambiguously separated;^ as was possible for the T- 
dependent measurements in alkanes^ and Gai^ However, 
this does not affect our main conclusion that the surface 
roughness, as probed by the DS, is entirely due to ther- 
mally activated capillary waves. 

The local density profile of liquid In is compared with 
that of Ga in Figure E^b) . Profiles derived from the ex- 
tremal parameters of fits to R{qz) and from extremal val- 
ues for temperature, surface tension and resolution func- 
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FIG. 4: (a) Comparison of the local real space density profile 
for liquid Indium ( — ) with the thermally averaged density 

profile for the same metal ( ). The averaged density profile 

is directly accessible by experiment and has been measured at 
170° C. (b) Real space local density profile for liquid Gallium 

( ) and liquid Indium ( — ) after the removal of thermal 

broadening. Densities p{z) are normalized to the bulk den- 
sities Poo of liquid Gallium and Indium. Inset: Decay of the 
amplitude of the surface-normal density profile for liquid Ga 
(open circles) and liquid In (filled circles). The lines repre- 
sent the fit of this decay of the surface layering to the form 
exp{—z/l) for liquid Ga ( — ) and In ( ). 



tion are indistinguishable on the scale of the figure. The 
amplitude of the first density peak for In is comparable to 
that of Ga. This reflects the fact that in the region near- 
est the surface, the density oscillations are determined 
primarily by ao, which has similar values for the two met- 
als. Further towards the bulk, layering is seen to decay 
faster for In than for Ga. We quantify this observation 
by the excellent fit of the peak amplitude of the electron 
density profile to the form exp{— z / 1) . The decay lengths 
I = 3.5A±0.6 for In and I = 5.5A±0.4 for Ga differ by 
an amount well outside of the experimental error. This 
is illustrated in the inset of Figure QJb), a semilogarith- 
mic plot of the maxima of the density oscillations as a 
function of distance from the first surface layer. The lines 
represent the fit of the peak amplitude to the exponential 
form and the slope represents the decay length. 

It is also instructive to compare the surface structure 
factor measured in the refiectivity experiments to the 
bulk structure factor measured by standard x-ray or neu- 
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FIG. 5: (a) Bulk (A) and surface (•) structure factors for liq- 
uid In. The bulk structure factor data are taken from Orton 
et al.^^ The solid line is a fit of the model explained in the 
text to the experimentally determined surface structure fac- 
tor. The broken line is a guide for the eye. The inset compares 
the decay of the longitudinal surface density oscillations (filled 
circles; solid line represents fit to an exponential) to the decay 
of the bulk pair correlation function (open triangles; broken 
line represents fit to exponential). The coordinate z of the 
X-axis represents the distance from the surface in the case of 
the decaying surface layering and the radial distance from a 
reference atom in the bulk liquid in the case of the decaying 
bulk pair correlation. The data for the bulk pair correlation 
function are taken from Waseda^ (b) Bulk (A) and surface 
(•) structure factors for liquid Ga. The bulk structure factor 
data are taken from Narten et al^ Solid lines are fits to re- 
fiectivity data, broken lines are guides for the eye. The inset 
compares the decay of the longitudinal surface density oscil- 
lations (filled circles; solid line represents fit to exponential) 
to the decay of the bulk pair correlation function (open trian- 
gles; broken line represents fit to exponential). The data for 
the bulk pair correlation function are taken from Waseda4- 



tron diffraction (Figure El^a): Indium; Figure El^b): Gal- 
lium) . The widths of the peaks shown are inversely pro- 
portional to the decay lengths of the corresponding cor- 
relation functions: the pair correlation function g(^r)^ 
of the bulk and the layered density profile p(z) of the 
surface. Plotting the amplitudes of g(r) and p{z) on a 
semi- log scale (inset of Figure EJ shows that they lie on 
straight lines. This indicates that the correlations de- 
cay exponentially, and the decay lengths are obtained as 
the negative inverses of the slopes. The inset plots show 
clearly that while the surface and bulk decay lengths are 
about the same for In, in Gallium the surface induced 
layering decays much more slowly than the bulk pair cor- 
relation function. This is consistent with the NFE nature 
of In. For Ga, by contrast, the tendency towards cova- 
lency may disrupt the hard sphere packing and enable 
surface correlations that differ from those in the bulk. 

Another piece of evidence suggesting that directional 
bonding in the melt might have an effect on surface in- 
duced layering stems from the ratio between the surface 
layer spacing and the bulk nearest neighbor distance. 
This nearest neighbor distance is taken from the anal- 
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ysis of the pair correlation function which is the Fourier 
transform of the structure factor depicted in Figure \^ 
Due to truncation problems and related subtleties with 
the data analysis, the bulk nearest neighbor distance can- 
not be directly inferred from the maximum in the bulk 
liquid structure factor and we use the nearest neighbor 
distance given by lida and Guthrie.'^® For liquid In, this 
ratio of surface layer spacing to bulk nearest neighbor 
distance is 2.69A/3.14A = 0.86, close to the value of 
^2/3 « 0.82 for the hard sphere packing that would 
be expected for a NFE liquid metal. For liquid Ga, this 
ratio is larger, 2.56A/2.78A = 0.92. In this case, for- 
mation of directional bonds, eventually leading to Ga2 
dimers, may prevent close packing. 

Tomagnini et al. recently considered a possible rela- 
tionship between surface induced layering in LM and the 
stability of crystal facets at metal surfaces. Although 
premclting of crystalline interfaces is quite common in 
non-metallic crystals, most metals have at least one close 
packed face which does not premelt, remaining solid up to 
the melting temperature^^. Tomagnini et al. use Molec- 
ular Dynamics simulations to demonstrate that when the 
period of the surface induced layering at the liquid sur- 
face, which they assumed to be 27r/qo where qq is the 
position for the peak in the bulk liquid structure factor, 
is commensurate with the distance between lattice planes 
along the normal to the crystal facet, this particular facet 
is stabilized and resistant to premelting. It would be in- 
teresting to correlate the premelting properties of Ga and 
In crystals with our results for the period of the surface 
layering and in particular, to examine if the effect of di- 
rectional bonding that affects the surface structure in the 
liquid affects premelting of the solid as well. 

Our diffuse scattering measurements show that the 
roughness of the liquid In surface can be attributed en- 
tirely to thermally excited capillary waves. This indi- 
cates that there are no other detectable inhomogeneities 
on the surface, such as might be caused by microscopic 
oxide patches. This observation is in concert with in- 
dependent experiments on organic Langmuir monolayers 
on water. If the monolayer is homogeneouaMi^, the DS 
can be described by Eq.l^UJ) with no adjustable parame- 
ters and no excess scattering is observed — just as is the 
case for metallic liquid In. On the other hand, if the or- 
ganic monolayer is compressed beyond its elasticity limit 
and becomes inhomogeneous, Ea. ljlOl) no longer describes 
the experimentally determined DS, and excess scattering 
is observed that must originate from sources other than 
thermally activated capillary waves44 This demonstrates 
the viabihty of diffuse scattering as a tool for studying 
surface inhomogeneities. This is of particular interest for 
investigations on liquid alloys where — depending on the 
type of alloy — a rich surface behaviour is expected, rang- 
ing from concentration fluctuations for alloys displaying 
surface segregation to critical fluctuations for alloys with 
a critical consolute point4& 

Controlled oxidation of the liquid metal surfaces re- 
veals a further striking difference between Ga and In. 



Although exposure of liquid Ga to air or to a large 
amount of oxygen leads to the formation of a macroscop- 
ically thick and rough oxide film, exposure to controlled 
amounts of oxygen under UHV conditions produces a 
well defined uniform 5 A thick oxide film.^'' This film 
protects the underlying bulk phase to a certain extent 
from further corrosion, as further dosage at low oxygen 
pressures (< 2 x lO^"' Torr) was found to have no effect 
on the oxide thickness or coverage fraction. This pasivat- 
ing of the surface is rather unusual for solid metals, Al 
being prominent among the few examples. Liquid In, un- 
like Ga, was not found to form a passivating oxide film. 
Instead, macroscopic clumps of oxide formed which did 
not wet the clean In surface. This corrosion mechanism is 
rather common for solid metal surfaces, the oxidation of 
Fe being the best known example. It is not clear to what 
extent this fundamental difference in the mechanism for 
corrosion is related to differences in the structure of the 
liquid surface and to what extent to the chemical affinity 
for oxidation of the respective liquid metal. 



VI. SUMMARY 

We have measured the x-ray reflectivity and small an- 
gle off specular surface diffuse scattering from liquid In 
at 170°C. Our results can be quantitatively explained 
by the convolution of thermally excited surface waves 
and a temperature independent surface structure factor, 
corresponding to theoretically predicted surface layering. 
The absence of excess diffuse scattering beyond that due 
to thermally excited capillary waves demonstrates that 
the liquid-vapor interface is homogeneous in the surface- 
parallel direction. The intrinsic layering profile of liquid 
In, obtained by removing the capillary wave roughen- 
ing, is compared to that previously reported for liquid 
Ga. For Ga, surface layering persists farther into the 
bulk than is the case for In. This may be attributed to 
directional bonding in the Ga bond stabilizing surface in- 
duced layering over a larger distance than is the case for 
liquid In. Further evidence suggesting a correlation be- 
tween the degree of covalency in the melt and the surface 
structure stems from the observation that the compres- 
sion of the surface layer spacing relative to the bulk near- 
est neighbor distance is close to the behavior expected for 
an ideal hard sphere liquid for In but considerably dif- 
ferent for Ga. Controlled oxidation of liquid In results in 
the formation of a macroscopic rough oxide, unlike the 
passivating microscopic uniform oxide film that forms on 
liquid Ga. 
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